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We present analytical and numerical results that demonstrate the presence of anomalous entan¬ 
glement behavior on the Dirac Billiards. We investigate the statistical distribution of the character¬ 
istic entangled measures, focusing on the mean, on the variance and on the quantum interference 
terms. We show a quite distinct behavior of the Dirac Billiard compared with the non-relativist 
(Schrodinger) ones. Particularly, we show a very plausible Bell state and a sharp amplitude of 
quantum interference term on entangled electrons left from the Dirac Billiards. The results have 
remarkable relevance to the novel quantum dots build of materials like graphene or topological 
insulators. 
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The entanglement is one the most fundamental effect 
on quantum mechanics, with no classical analog^. Two 
or more particles are entangled if they support a nonlocal 
correlation which cannot be acquired by the dynamics of 
classical mechanics^. Because its non-classical character¬ 
istics, the control of entangled states has attracted the 
interest of numerous science communities^. The tech¬ 
nological applications of the effect has a broad range as 
quantum computation, teleportation, telecommunication 

and cryptography^—. 

A large number of mechanism to entangle electronic 
particles, with or without interaction, can be found in 
the literatureAii, and the quantum chaotic devices are 
a promising option^. In a recent work, Beenakker 
et al . 9 proposed the possibility to entangle two non¬ 
interacting electrons using, as an orbital entangler, a 
quantum (Schrodinger) chaotic billiard or, as we will call 
henceforth, the Schrodinger Billiard (SB). They obtained 
the averages and variances of concurrence and entangle¬ 
ment. However the results were found to be nearly in¬ 
variant under the presence or absence of time-reversal 
symmetry (TRS). This means that the quantum interfer¬ 
ence corrections (weak-localization) of two arbitrary en¬ 
tanglement measures are approximately null, in contrast 
with another physical observables of electronic transport 
as conductance and shot-noise power. Nevertheless, the 
Ref. has argued that the two first moments of en¬ 
tanglement measures do not capture the full information 
about quantum dynamics, and the complete information 
about fundamental symmetries of nature emerges only on 
the distribution probabilities. More recently, the effects 
of tunneling barrier s 11 ! 12 on statistic of concurrence and 
joint probability distribution of concurrence and squared 
nor m —for SB were also studied. 

The novel materials, including Dirac materials 14 - 
(graphene and topological insulators), introduce new fun¬ 
damental symmetrie s 15 ! 16 that can affect the entangled 
electrons— In a recent investigation, the Ref.[0] an¬ 
alyzes how the sublattices or chiral symmetries^ (SLS) 
affect the electronic transport through a chaotic quantum 


(Dirac) billiar d 20 ! 21 , which henceforth we call chaotic 
Dirac Billiard (DB) 22 . Through this device, the wave 
functions of the electrons are described by massless Dirac 
equation of the corresponding relativistic quantum me¬ 
chanics, instead of Schrodinger equation. The two cat¬ 
egories of billiards show blunt differences on the corre¬ 
sponding electronic transport statistics, which leads us to 
suspect that a subtle difference occurs also on the quan¬ 
tum entanglement statistical moments. 

In this work, we analyze the concurrence and entan¬ 
glement statistics of two non-interacting electrons, with 
coherent phase transport, inside a chaotic DB. Two leads, 
both with two open channels, connects the DB with elec¬ 
tronic source and drain at the Dirac poin t 16 ’ 18 . We cal¬ 
culated the exacts distributions of concurrence (C) for 
systems with or without TRS. We show a drastic dis¬ 
similarity on the characteristic distributions of the DB 
compared with the SB results*— We highlight the two 
main difference between DB and SB focusing on plausible 
experiments. Particularly, first we show that it is more 
probable the production of maximally entangled elec¬ 
tronic state (C = 1 or Bell state) than of separable state 
(C. = 0) for experiment using a DB, precisely the opposite 
to the behavior of SB, as indicated in Fig. ©• Secondly, 
we show a weak localization correction of concurrence 
with the same order of magnitude of the mean (semi- 
classical) term, while for a SB it is approximately null^. 
Clearly, we indicate that the average and the variance of 
concurrence carry relevant information about quantum 
mechanics statistics of the DB, unlike what happens in 
SB. Our analytical results are in accordance with numer¬ 
ical simulation from the random matrix theor y — i—'— 

Scattering Model - The setup consists of a DB con¬ 
nected to two leads, both with two open channels. We 
consider the system in the absence of electron-electron 
interactions and at zero temperature. The DB is repre¬ 
sented by the massless Dirac Hamiltonian with SLS— 
The Hamiltonian satisfy the following anti-commutation 
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relationi^^ 


"H = -oyHoy, a z 


1 m 0 

0 — 1 m 


(1) 


The 'H-matrix has dimension 2 M x 2M, with 1 m denot¬ 
ing a M x M identity matrix. We can interpret the M 
number of l’s and — l’s in <j z as the number of atoms 
in the sublattices A and B respectively^, in a total of 
2 M atoms in the DB. The Hamiltonian model for the 
scattering matrix process can be written 

5(e) = 1 - 2 7 nW t (e-n + iTrVVW 1 ')" 1 W. (2) 


The 5-matrix has dimension 4x4, indicating a total of 
two open channels in each terminal, each one originated 
from A or B sub-lattices. The 2 M x 4 matrix W rep¬ 
resents all deterministic couplings of the DB resonances 
to the open channels of the two terminals. The scatter¬ 
ing matrix is unitary 5^5 = 1 due to the conservation 
of electronic charge. It is convenient to represent the 5- 
matrix as a function of transmission, t , and reflection, r, 
blocks as 


where t, t', r and r’ have dimension 2x2. From Eqs. © 
and ©, the 5-matrix also satisfies the relation 


5 — S Z 5^S Z , 


1 2 0 
0 - 1 2 ’ 


(4) 


at the Dirac point, zero energy (e = 0). Substituting the 
Eq.© in Eq. (fUjl, we conclude that r = r ', r 1 = r'l and 
t' = —if 

We follow the previously mentioned RMT model for 
the H matrix and, consequently, for the 5 matrix to 
investigate the entanglement between two no interact¬ 
ing electrons as described in Refs. [30]. Firstly, the con¬ 
currence of two electrons, after they left the mesoscopic 
device with SLS, can be written in terms of the trans¬ 
mission, tt\ eigenvalues, t\ and r 2 , both encoded in the 
scattering matrix, Eqs. © and © , as 

c = 2 Vn( 1-T,)r 2 (l —r.) 

n + t 2 — 2 tit 2 

The electronic states, after they left the DB, are sepa¬ 
rable or non entangled if C = 0, while, if C = 1, the parti¬ 
cles are maximally entangled (Bell state). For intermedi¬ 
ate values of C between 0 and 1, the states are known as 
non-separable or partly entangled^!. The entanglement 
and concurrence are related through the equation^ 


e(C) = h 


'l + y/T^C 2 ' 


( 6 ) 



FIG. 1: The distributions of concurrence ( 'Pp(C )) for the DB with 
(/3 = 1) or without (= 2) TRS. The Eqs. dill ) and dl2[ are depicted 
together with the corresponding ones for the SBi£. 


Probability Distributions - To obtain the distribution of 
DB’s concurrence, we start using the joint distribution, 
defined for the two characteristic eigenvalues, which was 
calculated in Ref. [23} for /? = {1,2}. The result can be 
written as 

Vpfai,<h) = sin(^>i T <j 6 2 ) sin(()>i 0 2 )|^ 

x sin^ _ 1 ( 20 i)sin^ _ 1 ( 202 ). (8) 


The relation between the variables fa and the transmis¬ 
sion eigenvalues ly is 7y = sm 2 {2<pi ), with fa € [o, 7 t/2], 
while normalization constants assume the values C\ = 1 
and c 2 = 6. 

The Eqs.© and © can be used to calculate the con¬ 
currence distribution Vp{C) following the definition 




C -2 


a/ti(1 - ri)r 2 (l - t 2 ) 
Ti + t 2 - 2tit 2 


(9) 


with (...) denoting the ensemble average performed with 
the eigenvalues distributions of Eq. ©. The Eq. © can 
be rewritten in the integral form as 


VpiC) = 


C — 2- 


tan(2<(>i) tan(2<^ 2 ) 


tan 2 (2</>i) + tan 2 (2^> 2 )_ 

x 'Pp(fa,fa)dfadfa (10) 


with 


h(x) = -x\og 2 (x) - (1 - x) log 2 (1 - x). 


The double integral over variables fa and fa in Eq. 
(7) (fTUl) can be performed through the transformation of 































FIG. 2: The distributions of concurrence ('P/ 3 (C)) and entanglement {Qp{e)) from Eas Jllt . d 1 21) and d 1 61) are depicted (continuous lines) for 
systems with (/3 = 1) and without (/3 = 2) TRS. The histogram is the numeric simulations of the Hamiltonian model with 10 6 realizations of the 
corresponding random »S-matrix, the Eq. 0. The analytical and numerical results agrees nicely. 


variables Zi = tan 2 (2 fa) and an appropriate expansion 
the delta function^ 


C-2 


tan(2</>i) tan(2 <j> 2 ) 
tan 2 (20i) + tan 2 ( 2 </> 2 ) _ 


with 


?±(C) = 


f±(C) = 


The Fig.JD shows the distributions of concurrence, 
Eqs. CD and (HU). for the DB. We plot in the same 
place the distributions for the SB, which were obtained 
= F+(C)z 2 6 [ Zl - /+(C)* 2 ]in Ref.®, for a direct comparison. The systems exhibit 

a quite distinct behavior. We observe, as the main char- 
+ J r -{C)z 2 S [zi — /-(C)22]acteristic of DB, a major probability to find the elec¬ 
trons in the maximally entangled states (C = 1), con¬ 
trary to the behavior of the SB. However, the preserved 
TRS supports electrons with the probability of separable 
states (C = 0) significantly increased, as depicted in the 
Fig. CD- The mean value of concurrence can be obtained 
from Eqs. CD and CD and renders 


(/±(c) + i) 2 y7^(c) 

±/±(C)T 1 
2-C 2 ±2 v / l-C 2 
C 2 ' 


After a some algebra, we obtain the following expres¬ 
sions for the probability distribution of concurrence with 
preserved TRI (/3 = 1) 


Vi (C) 


11 ( y/2 I VT—& 

21 C 2 | 2 V + 1-C 2 


( 11 ) 


x 


arccoth 



and, for broken TRI (3 = 2), it renders a simple expres¬ 
sion 


f 0.4675 /3 = 1 

\ 0.5708 J3 = 2 


(13) 


Using Ea. dTdl) . the weak localization of concurrence can 
be calculated with the same algebra and renders 


(C) wl = (C) 1 - (C) 2 « —0.1033. (14) 

Notice the weak localization term has the same mag¬ 
nitude order of the main (semiclassical) term, Eg. (1131) . 
while for the SB it is approximately null^. The TRI, 
unlike what happens in SB, has a strong influence over 
average of concurrence in DB. For the variance of con¬ 
currence we obtain 


W) 


Vi -c 2 

( 1 +C)( 1 -C 2 )' 


( 12 ) 


var [C\ 


0.1117 /3 = 1 

0.1034 (3 = 2 


(15) 






































































































4 


Like weak localization, the variance of concurrence lias 
the same order of magnitude of the its mean, Eq. tfldl) . 
Lastly, we obtain the distribution of entanglement Qp 
from Eqs. m and m using the following appropriate 
change of variables^ 


Q/J(e) 


hi(2) yi - C(e) 1 2 
C{e) arctanhf-^/l — C(e) 2 4 5 ] 


Me). 


(16) 


The analytical results for the distributions of entangle¬ 
ment are cumbersome, but it is depicted in the Fig. ([2]). 
The average of entanglement renders the results 


f 0.382 0 = 1 
\ 0.485 0 = 2 


(17) 


where 2 M is the number of resonances inside the DB, 
including both the two sub-lattices degrees of freedom. 
Also, A = 2MA/tt is the variance, related to the 
electronic single-particle level spacing, A. The W = 
(Wi,W 2 ) matrix is a 2 M x 4 deterministic matrix, de¬ 
scribing the coupling of the resonances states of the 
chaotic DB with the propagating modes in the two ter¬ 
minals. This deterministic matrix satisfies non-direct 
process, i.e., the orthogonality condition = -^S P ,q 

holds. Accordingly, we consider the relation cr z WE z = 
W, indicating the scattering matrix is symmetric (JU). We 
consider the system on the Dirac point, e = 0, and, to 
ensure the chaotic regime and consequently the univer¬ 
sality of observables, the number of resonances inside DB 
is large (.M ^>4). 


while its weak localization amplitude term is 


<£)„! = (e)i - <e> 2 « -0-10. 


For the variance of entanglement we obtain 


1 0.121 0 = 1 

{ 0.122 0 = 2 


(18) 


(19) 


Numeric Simulation - In order to confirm the analyti¬ 
cal results, the Eqs. ED and ED. we employ a numerical 
simulation using the 5-matrix formulation of Eq. (f2l) . The 
ensemble of Hamiltoninans satisfies the SLS of Eq. ([TD . 
The numerical simulation was developed in Ref. 0 . The 
anti-commutation relation, Eq. CD, implies a Hamiltonian 
member of the ensemble rewritten asi^ 

(20) 

Here, the T-matrix block of the H- matrix has dimension 
M x M. The Random Matrix Theory establishes that 
the entries of T -matrix can be chosen as a member of a 
Gaussian distribution^ 

P(T) cx exp |-|^Tr(r t T)|, (21) 


The numerical simulations produce the Fig.@, which 
shows the distributions of concurrence and entanglement 
obtained through 10 6 realizations compared with the an¬ 
alytical results, Eqs. ED. ED and ED. for systems with 
or without TRS. We use the T matrices, with dimension 
100 x 100 (M = 100), and the corresponding H matrices 
with dimension 200 x 200 (200 resonances). 

Conclusions - We present a complete statistical study 
of the entangled electronic measures. We show a very 
peculiar quantum behavior if the production of orbital 
entanglement is raised on the Dirac billiards quantum 
dots. The analytical expressions for the distributions of 
concurrence were obtained in the presence or absence of 
TRS. We compare the results for the Dirac billiards with 
previous ones of the Schrodinger billiards, Refs. 10, 
and found significant statistics differences between these 
devices. Clearly, we indicate that the average, the vari¬ 
ance of concurrence and entanglement carry relevant in¬ 
formation about quantum mechanics statistics of the DB, 
unlike what happens in SB. Finally, our analytical results 
are confirmed by numeric simulation from random matrix 
theory. 

This work was partially supported by CNPq, CAPES 
and FACEPE (Brazilian Agencies). 
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